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Abstract. A great deal of our understanding of star formation in the local universe 
> i has been built upon an extensive foundation of Ha observational studies. However, 

recent work in the ultraviolet (UV) with GALEX has shown that star formation rates 
(SFRs) inferred from Ha in galactic environments characterized by low stellar and gas 
densities tend to be less than those based on the UV luminosity. The origin of the 
discrepancy is actively debated because one possible explanation is that the stellar initial 
mass function is systematically deficient in high mass stars in such environments. In 
this contribution, we summarize our work on this topic using a dwarf galaxy dominated 
sample of ~300 late-type galaxies in the 1 1 Mpc Local Volume. The sample allows 
us to examine the discrepancy between Ha and UV SFRs using a statistical number of 
galaxies with activities less than 0.1 M yr _1 . A range of potential causes for such an 
effect are reviewed. We find that while the IMF hypothesis is not inconsistent with our 
observations, alternate explanations remain that must be investigated further before a 
final conclusion can be drawn. 
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1. The Original Motivation for this Study: Probing Starbursts in Dwarf Galaxies 

Since the late 70's and early 80's, a time when a number of foundational papers on 
the measurement of star for mation in external galaxies were published (e.g JKennicuttl 
19831 : iGallagher et alJll984h . the most common way of tracing star formation in nearby 



dwarf galaxies has been through the Ha emission line. And since that time, observa- 
tions of apparently isolated dwarf galaxies undergoing starburst activity have raised the 
question of whether the star formation histories of low-mass systems are more gener- 
ally characterized by burs t s rather than modes that are quiescent and c ontinuous (e.g., 
Gallagher & Huntedl 19841 : iHunter & ElmegreerJl200i iLee et alJl2009hh . 



An opportunity t o investigate this qu estion in a new way came with the launch of 
the GALEX satellite dMartin et alJl2005h . which enabled sensitive imaging of the far 
ultraviolet emission (FUV; ~1500 A) from hundreds of dwarf galaxies in the nearby 
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universe de Paz et al. 2007; iLee et al. I I2010I) . In combination with Ha measure- 
ments for a complete sample of galaxies, constraints on the average duration, frequency 
and amplitude of starburst episodes in the population can be computed, since the star 
formation timescales probed by the two tracers flank their expected duration (i.e., on 
the order of a dynamical time, roughly 10 8 yrs). To review, Ha nebular emission arises 
from the recombination of gas ionized by the most massive O- and early-type B-stars 
(M* > \1Mq). It therefore traces star formation over the lifetimes of these stars, which 
is on the order of a few million years. In contrast, the UV flux primarily originates from 
the photospheres of a wider range of O- through later-type B-stars (M* > 3M Q ), and 
thus measures star formation averaged over a longer ~ 10 8 yr timescale. 

Qualitatively, galaxy populations with predominantly bursty (or otherwise discon- 
tinuous) star formation histories should imprint a clear signature on the distribution of 
Ha-to-FUV flux ratios. Relative to populations where the activity is generally continu- 
ous, the average Ha-to-FUV ratio should be depressed. When the SFR is constant, the 
ratio assumes an equilibrium value when the birth of stars responsible for the FUV and 
Ha emission balance their deaths. Variations in the SFR over timescales on the order 
of ~ 100 Myr disrupt this equilibrium. In the time following a burst of star formation, a 
deficiency of ionizing stars develops as they expire relative to lower-mass, longer-lived 
B-stars that produce si gnificant amounts of UV emission, and the Ha-to-FUV flux ra- 
tio th us is lower (e.g., ISullivan et al.ll200(jl : ISullivan et al.ll2004l : llglesias-Paramo et aL 



2004). For predominantly bursty populations, determinations of the SF R based on stan 



dard linear conversions between the luminosity and the SFR (e.g., iKennicuttl 11998 



hereafter K98) should also find that the Ha-based SFRs tend be lower than those com- 
puted from the FUV, since the conversions are derived under the assumption of constant 
star formation. 

Previous UV studies reported hints of such a trend for low luminosity galaxies, 
and tentatively attributed it to a systematic increase in the prevalence of bursts in the 
recent SFHs of dwarf galaxies dBell & Kennicuttll200ll) . However, deviations from the 
expected ratio as a function of luminosity began to appear only when integrated SFRs 
less than about 0.1 M Q yr" 1 were probed (i.e., a few times lower than the SFR of the 
Large Magellanic Cloud), and just handfuls of galaxies with such low activities were 
included in prior work. Thus, our original objectives were to examine whether the 
trend persisted with an unbiased, statistical sample of star forming dwarf galaxies, and 
if so, to use the average Ha-to-FUV ratios to produce new constraints on the temporal 
variability of the integrated activity. 

Analysis with our dataset indeed revealed that the dwarf galaxies had systemat- 
ically lower Ha-to-FUV values than more luminous, higher surface brightness spiral 
galaxies. This seemed to confirm earlier suggestions of generally bursty star forma- 
tion in low mass systems, and standard population synthesis modeling showed that 
bursts which elevate the SFR of galaxies by a factor of >100 for a 100 Myr duration 
would reproduce the observed ratios. However, such large amplitudes appeared to be 
in conflict with other independent determinations of the recent star formation history 
(SFH), which showed that more modest enhancements (factor of ~5) were typical. We 
therefore were led to consider a wider range of potential causes for the low observed 
Ha-to-FUV r atios in dwarf gala xies. Here, we distill the main points of our analysis, 
as reported in Lee et al.l (12009 ah . Discussion of this topic also appears in a number of 
other contributions in these Proceedings, by Meurer; Boselli et al.; Eldridge; Johnson; 
Calzetti; and Pflamm-Altenburg; among others. 
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2. Data 



Our study is based upon data collected by the GALEX 11 HUGS (11 Mpc Ha UV 
Galaxy Survey) and Spitzer LVL (Local Volume Legacy) programs. The sample is 
dominated by dwarf galaxies, and is thus ideal for studying the nature of systems 
wi th low SFRs. I ntegrated Ha , UV, and mid- to far-IR flux cat alogs are published 



Kennicutt et a l. (2008), and lLee et ail d2010h . iDale et all (l2009h . respectively. De- 



tively complete to M B < -15 and M H i > 2 x 10 8 M Q at the edge of the 1 1 Mpc volume 



in 

tails on the sample selection, observations, photometry, and general properties of the 
sample are provided in those papers. A brief summary of the dataset is given here. 

Our total Local Volume samp le contains 436 object s. Galaxies are compiled from 
existing catalogs (as described in iKennicutt et all 120081) . and the selection is divided 
into two components. The primary component of the sample aims to be as complete as 
possible in its inclusion of known nearby star-forming galaxies within given limits. It 
consists of spirals and irregulars brighter than B = 15 mag within 1 1 Mpc that avoid the 
Galactic plane (\b\ >20°). These bounds represent the ranges within which the original 
surveys that provided the bulk of our knowledge on the Local Volume galaxy popu- 
lation have been shown to be relatively complete, while still spanning a large enough 
volume to probe a representative cross section of star formation properties. The sec- 
ondary component of the sample consists of galaxies that are within 1 1 Mpc and have 
available Ha flux measurements, but fall outside one of the limits on brightness, Galac- 
tic latitude, or morphological type. It is a composite of targets that were either observed 
by our group as telescope time allowed, or had He fluxes published in the literature. 
Subsequent statistical tests, as functions of B-band apparent magnitudes and HI fluxes 
(compiled from the literature), show that the subset of galaxies with \b\ > 20° is rela- 

V compl 
dLee et alJl2009bh . 

Through a combination of new narrowband Ha+[NH] and i?-band imaging, and 
data compiled from the literature, integrated Ha fluxes are available for over 90% of 
the total sample. The new narrowband imaging obtained by our group was taken at 1-2 
m class telescopes in both hemispheres, and reached relatively deep point source flux 
and surface brightness limits of ~2 x 10~ 16 ergs cm -2 s -1 and ~4 x 10~ 18 ergs cm -2 
s -1 arcsec -2 , respectively. 

Subsequent GALEX UV imaging primarily targeted the \b\ > 30°, B < 15.5 subset 
of the sample. The more restrictive latitude limit was imposed to avoid excessive Galac- 
tic extinction and fields with bright foreground stars and/or high background levels for 
which observations would be prohibited due to GALEX's brightness safety restrictions. 
Deep, single orbit (~1500 sec) imaging was obt ained for each galaxy , following the 
strategy of the GALEX Nearby Galaxy Survey dGil de Paz et al.ll2007h . GALEX ob- 
servations for a significant fraction of the remaining galaxies beyond these limits were 
also taken by other GI programs. Overall, GALEX data are available for ~ 90% of the 
sample. 

Finally, Spitzer IRAC mid-infrared and MIPS far-infrared imaging was also ob- 
tained for the \b\ > 30°, B < 15.5 subset of the sample through the Local Volume 
Legacy program. The far-infrared photometry provide attenuation corrections for our 
analysis. 
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3. Results 

Our main results are shown in Figures Q] and |2] Figure Q] plots the FUV luminosity 
against the Ha luminosity, while Figure[2]shows the Ha-to-FUV flux ratio as a function 
of the Ha luminosity. A systematic decrease in the Ha-to-FUV ratio with decreasing 
luminosity is evident. At high SFRs (> 0.1 M G yr -1 ), the ratio is constant, and the value 
of SFR(Ha)/SFR(FUV) is ~0.75. Within the uncertainties, this value is consistent with 
a zero offset from the expected equilibrium value for constant star formation. However, 
by SFRs of about 0.003 M Q yr -1 , the average Ha-to-FUV flux ratio is lower by a factor 
of two, and at the lowest SFRs probed (10~ 4 M Q yr -1 ), the average deviation is about a 
factor of ten. 

Best-effort attenuation corrections have been applied to the values plotted in Fig- 
ures Q] and |2] (i.e., based on Balmer decrements and total infrared to FUV ratios when 
available, and on empirical scaling relationships otherwise). However, the trend is al- 
ready evident prior to these corrections; it thus should be robust to the effects of dust 
since the FUV should be more severely attenuated than the nebular Ha emission. 



4. Understanding the Systematic Decline in L(Ha)/L(FUV) 

Again, our analysis confirms previous indications of systematically lower Ha-to-FUV 
ratios in dwarf galaxies based on small datasets (~ 20 galaxies with SFRs < 0.1 M Q 
yr -1 ), and improves the sampling of this regime by an order of magnitude. Our re- 
sult is also corroborated by a number of concurrent, independent studies (though with 
fewer numbers of dwarfs), and the trend has als o been reported as a function of decreas- 
ing optical surface brightness and stellar mass dMeurer et al.ll2009l : iHunter et~aL 1 120101 : 



BosellietaL||2009). 



As discussed above, our initial assumption was that such a trend would be an 
indication of an increased frequency of starburst activity in the dwarf galaxy popula- 
tion, and that modeling the decline of the average ratio would yield new constraints 
on the characteristic durations, frequencies and amplitudes of the bursts. Examina- 
tion of a Starburst99 model grid spanning a range of burst parameters computed by 



Iglesias-Paramo et a l. (2004) showed that the factor of two offset observed at SFR-0.003 
M Q yr -1 could be reproduced by bursts with amplitudes of >100 lasting for 100 Myr. 
However, such large amplitude bursts appear to be in conflict with other observational 
constraints. The most direct constraints are provided by studies which reconstruct star 
formation histories from resolved observations of stellar populations in nearby low- 
mass systems (e.g. JWeisz et al.ll2008l : iMcOuinn et~ail2009l : iMcOuinn et aO2010h . The 
typical burst amplitudes found in such studies range from a few to ~10, an order of 
magnitude smaller than the factor of 100 bursts required by the models to depress the 
Ha-to-FUV flux ratios by a factor of two. The 11 HUGS sample itself also provides 
a statistical constraint on the average dwarf galaxy starburst amplitude via the ratio of 
fraction of star formation (as traced by Ha) concentrated in starbursts to their number 
fraction dLee et al.ll2009bl) . The high degree of statistical completeness of 11HUGS 
makes this calculation possible for galaxies with Mb < -15, and again, the burst am- 
plitude is found to be relatively modest (~4). 

We were thus led to consider a wider range of explanations for the trend, the most 
viable of which are: 
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Figure 1 . Comparison of FUV and Ha luminosities (corrected for internal dust 
attenuation), with axes indicating the corresponding Ha and FUV SFRs, com- 
puted from the linear conversion recipes of K98. The solid line represents a 
one-to-one correspondence between the SFRs. The top panel show the data with 
predictions from models which perfo rm population synthesis by random sam- 
pling of the IMF dTremonti et al J 120071) . The median-predicted values (solid line) 
are shown along with values at the 2.5 and 97.5 perce n tile p oints (dotted lines). 
Predictions from IGIMF mod el of iRroup a & Weidner ( 2003]), as computed by 
Pflamm-Altenburg et al. (2007) are shown with the data in the bottom panel. The 
bottom and left hand axes are shown in gray to signify that the K98 SFR scales 
would not be valid at low SFRs under the assumptions of the IGIMF model. 



(1) a leakage of ionizing photons from the low density environments which 
are characteristic of the galaxies where low Ha-to-FUV ratios are ob- 
served; 
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Figure 2. Analogous to Figure 1, but for the ratio of Ha-to-FUV 
SFRs/luminosities (corrected for internal dust attenuation) as a function of the Ha 
SFR/luminosity. The shaded band represents the range of Hc-to-FUV ratios pre- 
dicted by commonly used stellar population models for constant star formation. The 
dashed line shows a linear least squares fit to the data. Circled points represent galax- 
ies where the Ha flux may be underestimated because the narrowband imaging did 
not wholly enclose the galaxy. The top panel shows the deviations expected from 
random sampling of the IMF, while the bottom panels shows the predictions of the 
IGIMF model. 



(2) low probabilities of forming ionizing stars when the integrated SFR of 
a galaxy is low; 

(3) an IMF which is systematically deficient in the highest mass stars in 
low density environments. 

The first possibility is notoriously difficult to observationally constrain, though 
the few existing studies of Lyman continuum escape from present-day gala xies have 
all found upper limits <10% (e.g.. iLeitherer et ailll995l : iBergvall et al.ll2006h . These 
detection experiments have been performed on starbursting galaxies where leakage is 
thought to most likely occur. Although it is difficult to conclusively rule out this possi- 
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bility, given the reported upper-limits, it seems improbable that ionizing photon losses 
could approach the levels required to explain the factor of two discrepancies observed, 
particularly in the quiescent, gas rich dwarfs that dominate our sample. 

A related explanation is that low gas densities lead to extremely diffuse Ha emis- 
sion that falls bene ath the detection threshold of standard narrowband observations 
dHunter et alj|20ldh . In this case, the ionizing photons travel far from their parent HII 
regions, but do not to escape the galaxy entirely. To check whether this could be an issue 
for our data, we performed the following simple calculation. The potentially missing 
Ha luminosity is computed assuming that the Ha star formation rate (SFR) should be 
consistent with the SFR inferred from the UV continuum flux. The surface brightness 
of this emission is then estimated with the additional assumption that the nebular flux is 
uniformly distributed over the area of the UV disk. We found that, under these assump- 
tions, the missing Ha should have been unambiguously detected in our observations. 
However, if the nebular emission instead extends to twice the area of the UV disk (this 
generally encloses the periphery of the HI gas distribution, and thus serves as a lim- 
iting case), then the emission will fall below our sensitivity limits. This issue is one 
that requires further examination. To do this, we recently obtained deep narrowband 
(~8A) imaging of a few dwarf galaxies with the Maryland-Magellan Tunable Filter on 
the Inamori Magellan Areal Camera and Spectrograph (IMACS) at the Magellan 6.5m. 
The data will allow us to observe Ha down to ~4 xl(T 18 ergs s _1 cm^ 1 arcsec~ 2 — 
sensitive enough to detect the surface brightnesses expected if the "missing Ha" is in 
fact spread over the entire area of the gas disk. 

The second possibility is that random sampling of a universal IMF in the regime of 
ultra-low integrated SFRs, may lead to an apparent deficiency or absence of Ha, even 
when there is on-going star formation, because the probabilities of forming high-mass 
ionizing stars are low. We use Monte Carlo simulations to model the impact of such 
stochasticity on the Ha-to-FUV flux ratio, and compare the results to the data in the top 
panels of Figures Q] and [2 The median-predicted values (solid line) are shown along 
with values at the 2.5 and 97.5 percentile points (dotted lines). While stochasticity 
clearly does have an effect, its impact does not appear to be large enough to explain the 
observed trend by itself. 

Finally, we also considered the possibility that the IMF is systematically deficient 
in the highest mass stars in the low density environments of dwarf galaxies. In particu- 
lar, w e compare the model of the Integrated Galactic IMF (IGIMF; iKroupa & Weidner 
20031 : also see Pflamm-Altenburg; and Weidner, in these Proceedings) to the data in the 
bottom panels of Figures Q] and [2] The model assumes that the most massive star that 
can form depends on its parent cluster's mass in a deterministic (rather than a proba- 
bilistic) manner, while the mass of the most massive cluster that can form is dependent 
on a galaxy's integrated star formation rate. This model is able to reproduce the trend 
observed in our data. 



5. Discussion and Next Steps 

While the consistency between the IGIMF model and the data is intriguing, alternate 
explanations have not yet been explored to the point where they can be conclusively 
eliminated. Clearly, the true fate of the ionizing photons in low density gas will be 
difficult to determine, though we hope that our new ultra-deep tunable filter imaging 
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will provide new insight on this issue. There was general consensus at the meeting that 
this must be investigated further. 

Whether the starburst scenario, or some class of fluctuating star formation histo- 
ries, can be ruled out as the cause of the low Ha-to-FUV flux ratio was a question 
that spawned much greater debate (e.g., contributions by Boselli et al. and Meurer 
in these Proceedings). Our position has been that the modeling of star formation his- 
tories warrants further examination, particularly in conjunction with the modeling of 
stochastic sampling of the IMF. Although we currently find that the large burst ampli- 
tudes required to reproduce the trend are in conflict with other observational constraints, 
stochasticity may amplify the effects of bursty or non-uniform "flickering" SFHs on the 
Ha-to-FUV ratio, and reduce the required amplitudes. Fortunately, a new generation 
of stellar population synthesis models that incorporate random sampling of a standard 
IMF has been recently developed by several groups (e.g., see contributions by Eldridge; 
Fumagalli; and Cervino in these Proceedings). Such models are much better suited for 
probing the properties of systems where the IMF may not be fully sampled, compared 
to existing models such as Starburst99 which use fully -populated ~10 5 M Q quanta as 
the basis for the synthesis. Overall, it seems possible that some combination of all of 
the effects described above may conspire to produce the observed systematic. 

Although the origin of the trend is currently unclear, there is one reasonable con- 
clusion that can be drawn from our study: that FUV luminosities provides more accu- 
rate SFRs for individual galaxies with low total SFRs and low dust attenuations. The 
FUV emission should be less prone to stochastic effects from sparse sampling of the 
upper end of the IMF, and to possible uncertainties in the fate of ionizing photons. With 
the demise of the GALEX FUV detector however, there is no longer a facility capable 
of obtaining such data, and for the foreseeable future, new measurements of the SFR in 
dwarf galaxies will need to be based on Ha emission. The role of the new stochastic 
stellar population synthesis models is thus critical, as they will be necessary to inter- 
pret the most likely value and possible range of the SFR that a given Ha luminosity 
corresponds to. 
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